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The polarization transfer coefficient Ky for the 2 H(d^ jj) 3 H reaction at energies i?i a b. = 30, 50, 70, 
and 90keV is calculated using the four-body Faddeev-Yakubovsky equation. We employ a separable 
type of the Paris potential. The experimental data [Phys. Rev. C64, 04760f (200f )] agrees with our 
solution at the scattering angle ^c.m. = 0°, and there is no discrepancy pointed out in that paper. 
We found that the is sensitive to the incident deuteron energy at this angle. 
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Since Kulsrud et al. proposed the idea of a neutron- 
lean fusion reactor 0, low energy 2 H(d, p) 3 H and 
2 H((i, n) 3 He reactions have attracted much attention. 
Subjects of discussions are mainly about suppression of 
the reaction by the polarized beam and target deuterons 
in the same direction (d + d), to that by the unpolar- 
ized ones (d + d). Experimentally there is no direct data 
for the d + d reaction, however new data of the cross 
sections by the unpolarized 2 H(d,p) 3 H and 2 H(d, n) 3 He 
reactions J2J and the vector and tensor analyzing powers 
by the polarized beam || Q at several lOkeV energies. 

On the other hand there are many theoretical a 
proaches to the suppression problem, with models [B [~ 
and analyses of the experimental data ^ 0, |l0|]. But 
their solutions are all different each other depending the 
applied models or analyses methods, and there is no con- 
clusion. 

To obtain more accurate solution, we employ the four- 
body Faddeev-Yakubovsky (FY) equation |ll|] which is 
able to calculate the system in microscopic way. The FY 
theory is a direct extension of the three-body Faddeev 
theory to the four-body system, and it includes all kine- 
matics of four nucleons. There are many vigorous studies 
in this field. For recent example, the analyzing powers 
for 4N system Jl2l [l3| O), and the cross section and low 
energy resonance for Ti^He elastic scattering jl4|, are 
reported. And the FY equation is a ppli ed to the 4 He 



bound state |16|] and hypernuclei [y_7|, |18| . In regard to 
the 4 He ground state, our numerical result for the binding 
energy agrees with that by Fonseca |l9) within a limited 
number of state channels and employing Yamaguchi po- 
tential. 

In this decade, our investigations were mainly on the 
analyses of the experimental data of the cross section 
and analyzing powers, together with prediction of the 
suppression ratio j20], ||, |3[ ||, ||, |(|. It began from 
Ref. [M and estimation of the suppression by the invari- 



ant amplitude method 21 is compared with our predic- 
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tion using the FY equation in Ref. p2[ . The final result 
employing the Yamaguchi-type nucleon-nucleon (NN) in- 
teraction are presented in Ref. |^3|, |24j . And the first step 
using the realistic potential is discussed in Ref. J25| . We 
think that a study of the S-factor is one of our important 
subject 1 26 to present a data for the astrophysics. 

Recently Katabuchi et al. measured the polarization 
transfer coefficient K^' for the 2 H(d,p) 3 H reaction [p7| . 
The incident deuteron laboratory energy (Ed) is 90keV, 
however, the data included those for all energies between 
- 90keV with the energy loss in the target. The accep- 
tance angle of the detector was ±6.5° around the scat- 
tering angle 0° in the laboratory system. The obtained 
data was K$ = 0.09 ±0.10. 

In the same paper, they reported numerical results of 
the Kjj using the amplitudes determined by Lemaitre 
and Schieck S , where the amplitudes are expressed as a 
production of the Coulomb penetrability function and 
"internal" reaction amplitudes in the NN interaction 
without the Coulomb force. The "internal" amplitudes 
are treated as parameters independent of the incident 
deuteron energy, and fitted to many experimental data 
for the d + d reactions. The energy dependence is kept 
merely in the Coulomb penetrability function. Obtained 
Ry has very small energy dependence in all scattering 
angles, and the mean value over £"i a b. is different from the 
measured data by 1.5 standard deviation. Thus they con- 
cluded that the other models should be compared with 
the data. 

Complying to the discussion, we calculate the at 
energies i?i a b. =30, 50, 70, and 90keV using the FY equa- 
tion with following conditions. As the NN interaction, 
the Ernst, Shakin, and Thaler separable expansion |2|] 
of the Paris (PEST) potential |Q are employed for the 
following two-body states: % (3), 3 Si- 3 Di (4), 3 P (2), 
(2), 3 Pi (2), and 3 P 2 (3). The numbers in the braces 
are ranks for each states. In the [3+1] sub-amplitude, the 
three-body total spin and parity j 7r considered are within 
the range 1/2* to 7/2*. For the [2+2] sub-amplitude, all 
possible quantum states, which are automatically driven 
by two-body states, are taken into account. As for the 
four-body total spin and parity states J* these are taken 
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from 0* to 4^ and the convergence is confirmed. 

The [3+1] and [2+2] subamplitudes are expanded in 
separable form adopting the energy dependent pole ex- 
pansion (EDPE) |30| method with following modifica- 
tion. It is convenient to describe an outline of the tech- 
nique with the expressions in Ref. pl[ . The original 
EDPE method requires to solve the eigenequations (15) 
with fixed energies E x and E y which are usually chosen 
below the d + p + n threshold Energy. We adopt the nu- 
merical results of the binding energies for the triton as E x 
and two deuterons as E y in our calculation, and they are 
used in all states of the [3+1] and [2+2] subsystem. The 
form factors and propagators are obtained from the third 
through sixth equations in Eqs. (14) with required energy 
z on the FY calculation. In case for the mixed channels 
of the repulsive and attractive forces or the multi-rank 
separable potentials in the 2-body subsystem, Eqs. (15) 
have many complex solutions and the convergence of the 
FY calculation become poor. The modification intends 
to improve this difficulty. 

First the 2-body propagator r is diagonalized with a 
unitary matrix U as 



UtU~ 



(1) 



Next all the positive diagonal elements are set by negative 
but the same absolute values as 
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Thus the reduced propagator t r is obtained with U as 



(3) 



and t r is used in the eigenequations (15) except for r. In 
this case all the eigenvalues and eigenfunctions become 
real, and the FY calculation keeps good convergence. Af- 
ter obtaining the form factor, the [3+1] and [2+2] prop- 
agators are calculated with the original r. Details of this 
technique will be described on the other paper. 

As for the Coulomb force, exact calculation requires 
to distinguish the proton and neutron, and to include 
the pure Coulomb potential between two protons. How- 
ever nobody is successful to carry out this method in 
the FY integral equation within the momentum repre- 
sentation for the difficulty of the logarithmic singularity. 
An approximation method is proposed in Ref. g but we 
adopt a different way of which basic idea is described 
as "Traditional" in Ref. p6| . In this approximation, the 
FY equation for 4-identical particles is solved without 
Coulomb force, and obtained amplitudes are regarded 
as solutions for a mixed channel 2-body problem where 
the two deuterons in the initial channel and proton and 
triton in the final are treated as point particles. The 
scattering amplitudes are described the operator of the 
T-matrix put between the plain wave functions for the 
initial and final channels. However in case for existence 
of the Coulomb repulsion, the Coulomb functions ought 
to be applied except for the plain waves. We modify this 



approximately with multiplying the amplitude by simple 
factors as 
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where Mpji is a reaction amplitude without Coulomb 
whose t and I' are the angular momentum in the ini- 
tial and final channel, respectively. PoS and (5^s are the 
damping effects and phase shifts corrections, respectively, 
where without and with prime indicates the initial and 
final channel, respectively. These factors are expressed as 
the ratio of the Coulomb wave function to the plain wave 
which is applied the limiting value of the radial valuable 
r to be since the range of the NN interaction is much 
smaller than that of the Coulomb force. 
Thus we obtain Pq as 
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where ip c ( r ) an d 4>k( r ) ar< 3 regular functions of the 
Coulomb and plain waves, respectively, r) — 0.1581/ y/Ei 
is the Sommerfeld parameter where Ei is a center of mass 
(cm.) energy in the initial channel. Pq has the same ex- 
pression except for changing r\ into rf — 0.1369/ \[E~f 
where Ef is a cm. energy in the final channel. It is 
easy to obtain the Gamow factor from Pq in case for 
small energy which is described for instance in Ref. |52| . 
This factor is cancelled in due to independence of 
the partial waves. 

The Coulomb phase shift 8g is defined in the same man- 
ner as 



argT(£ + 1 + irj) — arg 



lim 



+o cj) e (r) 
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where ipe(r) and 4>e(r) are regular functions of the 
Coulomb and plain waves, respectively, represented by 
the partial wave decomposition. Si is represented in the 
same expression except for changing i and r\ into £' and 
rf, respectively. It is first introduced in Ref. |33| for a 
two-body system and applied it to the three-body sys- 
tem in Ref. @. 

One may think that the modification will be improved 
by changes of the functions V'c(r) and (^fc(r) in Eq. (||) for 
ijje(r) and <pe(r). However in this case, the P-wave com- 
ponents in the initial channel are too much enhanced and 
the results for the differential cross sections and analyz- 
ing powers become worse. Then we did not adopt it. 

FIG. shows our solution of together with the 
data in Ref. |27j . They are almost covering the error bar 
at #c.m. = 0°, thus the experimental data is described 
with the FY equation. The calculations are numerically 
converged enough for this discussion. The difference is 
about 2% between the present results and approximated 
ones including only 1 So, 3 Si- 3 Di states in the two-body 
subsystem. We expect that the influence to the results is 
small with inclusion of the higher spin states. 

Once we reported similar calculations in Ref. (24|, p5[ , 
however, the graphs for the show our results without 
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FIG. 1: The polarization transfer coefficient Ky at I?i a b. = 
30keV (solid line), 50keV (dashed line), 70keV (dotted line), 
and 90keV(dash-dotted line). The measured data (solid cir- 
cle) in Ref. |27| is compared with our numerical results. The 
triangle dot snows the mean value over the incident energies 
of their calculational results in the same paper. 
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FIG. 2: Same as for FIG. |Tj but for neglecting the Coulomb 
modification. 



the Coulomb modification. Then the discussions in the 
papers were based on these wrong results. 

We found two discrepancies between our results and 
those in Ref. ^7j. 1) Ours are sensitive to the incident 
deuteron energy at around # c . m . = 0°, where the larger 
the energy becomes, the smaller the value is. Then we 
expect if the target is as thin as the beam going through, 
the measured data will be small. 2) There are peaks 
at around 6> c . m . = 60°-90°. We suggest to measure at 
around those angles. 

To find out an origin of the sensitivity to the en- 
ergy at 6 c . m . = 0°, we compare the results for without 
Coulomb modification (FIG. ||) with the original (FIG. 
|l|). The former case is negative at around 9 c . m . = 0° 
and positive at around c . m . = 180°, which is totally 
different from the latter. And the difference between 
Eub. = 30keV and 90keV for the former is smaller than 



FIG. 3: Same as for FIG. [j] but for disregard of the P-wave 
components in the initial channel. 
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FIG. 4: Same as for FIG. [l] but for neglecting both of the 
P-wave components in the initial channel and Coulomb mod- 
ification. 



the latter by about 1/3 at 6 c . m . = 0°. Thus we think that 
most of the energy dependence is due to the Coulomb 
force. This is against to the numerical results in Ref. p7[| 
since the energy variation of their amplitudes depends 
only on the Coulomb penetrability function. 

We make another check on a connection between the 
energy dependence and partial waves. Since S- and P- 
waves are the main components in the initial channel 
of low energy 2 H(d,p) 3 H reaction, becomes almost 
isotropic by disregard of the P-wave components from the 
reaction amplitudes (FIG||). In this case the Coulomb 
effects are canceled within our modification by multipli- 
cation of an amplitude and one's complex conjugation. 
Therefore the difference is small between FIG. || and |] 
which is neglecting both of the P-waves and Coulomb. 
Then we find with comparison FIG. |l| - 1| that the energy 
dependence in FIG. |l] comes from the Coulomb modifi- 
cation which affects to the mixing term of S- and P-wave 
components in the initical channel. 
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In case for neglecting P-waves, energy dependence at 
$c.m. = 0° is small and our solutions are close to that by 
Katabuch et al. using the amplitudes in ref. [||. It is 
interesting to us and further investigation is in progress. 

The author appreciates helpful discussions with profes- 
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